# -*- coding: utf-8 -*-

"""
@File    : 01_栅格文件展示.py
@Author  : fungis@163.com
@notice  :
"""
# 导入gdal,注意导入的名称
from matplotlib import pyplot as plt
from osgeo import gdal  # 或者直接用import gdal

raster_path = r'./data-use/tif/NDVI_201405_bm.tif'

# 读取影像
dataset = gdal.Open(raster_path)
# 打开波段1（注意:用索引1，而不是0，来获取第一个波段）
band = dataset.GetRasterBand(1)

band_array = band.ReadRaster()  # 一维像元矩阵
location_array = band.ReadAsArray()  # 二维行列像元矩阵
# cmap='gray'以灰度图形式展示
plt.imshow(band.ReadAsArray(), cmap='gray')
plt.show()
